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Abstract 

Fireballs are an essential part of any cosmological 7-ray burst. We derive here a 
stability criterion for fireballs and show that fireballs are Rayleigh- Taylor unstable in any 
region in which the entropy decreases outward. The instability begins to operate when 
the fireball becomes matter dominated. Among the possible implication of the instability 
are: (i) Conversion of a fraction of the radiation energy to a convective energy expressed 
in the motion of bubbles relative to each other, (ii) Penetration of fast bubbles through 
slower ones and creation of high 7 regimes which are essential for efficient conversion of 
the energy to 7-rays. (iii) Formation of rapid variation (of the scale of the bubbles) in the 
observed 7-rays. 
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1. Introduction: Fireballs and Grbs 

Regardless of the nature of their source, cosmological 7-ray bursts (GRBs) inevitably 
involve the creation of an expanding "fireball": an extremely optically thick plasma of 
electrons and positrons. Even if the energy is initially entirely in the form of gamma- 
radiation, electrons and positrons will be spontaneously created by the reaction 7 + 7 — > 
e + + e~ (Schmidt 1978, Cavallo & Rees 1978). The understanding of fireball evolution is 
therefore crucial to the understanding of the observed GRBs. 

The current picture of fireball evolution is based on smooth spherical models (Good- 
man, 1986, Paczyhski 1986, Shemi & Piran, 1990, Piran, Narayan & Shemi, 1993, Meszaros 
and Rees 1992, Rees & Meszaros, 1992) with a possible inclusion of shocks (Narayan, 
Paczyhski & Piran, 1992; Rees & Meszaros 1994). In this letter we show that fireballs 
may be Rayleigh- Taylor unstable. The instability develops during and after the transition 
from radiation dominated to matter dominated fireball and therefore causes the evolution 
to differ from that predicted by a smooth spherical model. The existence of the instabil- 
ity may ease the constrains on the baryonic rest mass allowed to be entrained within the 
fireball, and may also have direct observational implications. 

2. Fireball Evolution, Asymptotic Laws and Shocks Waves 

Piran, Shemi & Narayan (1993) have analyzed the evolution of a fireball in terms of the 
evolution of its individual shells. This descriptions holds for an impulsive fireball, in which 
all the energy is released instantaneously (Goodman, 1986), as well as for a quasi-steady 
radiative wind, in which the energy is released over a period longer than the dynamical time 
(Paczyhski, 1986). Under certain conditions, it is also valid for non spherical fragments of 
a spherical shell (Piran, 1994; Narayan & Piran, 1994). 

Each shell is characterized by a radiation energy density e r , a baryon number density n 
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and a relativistic Lorentz factor 7. The evolution may be divided into two phases. Initially 
the shell is radiation dominated, 77 = e r /m p nc 2 ^> 1, and the evolution is described by: 
e r oc .R -4 , n oc -R -3 , 7 oc R and 1] oc R~ x (R is the fireball radius). This description 
is valid after a rapid initial acceleration that gives the shell relativistic velocity. Since 
77 decreases the shell eventually becomes matter dominated. In the limit of 77 <C 1 the 
baryons coast with: 7 ~ 77^ (the index i denotes initial quantities) and n oc R~ 2 . The 
radiation continues to cool with 77 oc R~ 2 / 3 . The transition between radiation and matter 
dominated evolution occurs roughly at R = R v = Rirji, where 77 ~ 1. 77 and 7 are related 
by 7(477/3 + 1) — 77/37 = 771 + 1 throughout the shell evolution. 

The fireball becomes optically thin (at R = R T ) during the radiation dominated phase 
only if 77 is extremely large (typically > 10 5 ). In this case the released radiation is thermal 
with T Ti{R/Ri) (Due to the relativistic expansion an observer at rest will see blue 
shifted radiation with T obs m 7T m Tj). If 77^ is substantially smaller the shell becomes 
matter dominated (77 ~ 1) before it becomes optically thin, and only an insignificant 
amount of radiation energy is left at R ^ R T . However, not everything is lost. At a much 
larger radius the baryons may interact with the ISM to produce a SNR like shock (Meszaros 
& Rees 1993; Meszaros, Laguna, & Rees, 1993) The kinetic energy is then converted back 
to radiation at the shock and it produces a GRB if 7>10 3 . The last condition poses a 
strong limit on the baryonic load: M < 10~ 6 M Q (£'/10 51 ergs), which might be difficult to 
satisfy. This stringent limit is the most serious open question that confronts cosmological 
GRB models today. 

This simple description breaks down when the model predicts shell crossing. Consider 
two adjoint shells denoted h and / with Rh < Ri separated initially by a distance Ar. When 
the shells become matter dominated they coast at 7^ ~ 77^ and 7; ~ rji t . If 77^ < 77^ the 
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inner shell moves faster and overtakes the outer one at R m + AR SC « R m + r] 2 Ar. AR SC is 
almost independent of rjh (provided that rjh 3> rji). Shell crossing signifies the breakdown 
of the simple description given in the previous paragraphs (in fact it breaks down slightly 
earlier) and the formation of shock waves. In general, the expansion is smooth if drji/dr > 
and shocks appear if drji/dr < 0. Narayan et. al. (1992) and recently Rees & Meszaros 
(1994) suggested that such shocks could convert a fraction of the kinetic energy back to 
radiation and, if this takes place at R > R T , produce GRBs even with a modest r\ fireballs. 

3. Stability Analysis 
3.1 Short Wavelength Limit and the Isobaric Approximation 

Large scale deviations from spherical symmetry lead, as mentioned in §2, to quasi- 
spherical local expansion. We shall therefore be interested only in the short wave-length 
limit in which the transversial linear scale of the perturbations, L, is small compared with 
the fireball's thickness -R/7. 

The baryon's thermal energy is always negligible in a fireball. Hence the pressure is 
radiation dominated, p = aT 4 /3, and the proper energy density is: e = nmc 2 + aT 4 . The 
speed of sound, c 2 s = c 2 /3(l + 3/4?y), is therefore comparable to the speed of light as long 
as r\ is not much smaller than unity. The time required to restore pressure equilibrium 
in a perturbed fluid element of size L (in the fluid's rest frame) is ~ L/c s ~ L/c. In the 
observer frame it equals 7-L/c. For L < R/7 this is shorter than the time for changes in 
the unperturbed solution, ~ R/c. Thus, for L < R/7 the perturbations are isobaric, i.e. 
we may assume that the pressure perturbation vanishes. (This assumption holds provided 
that the typical growth rate of the instability is not much faster than the typical rate of 
changes in the unperturbed solution). 
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3.2 Displacement Evolution 



We analyze the time evolution of a small fluid element displaced from its original 
trajectory, ro(t). Quantities along ro(t) are denoted by subscripts (70 (£), Po(^), eo(£) 
etc.), while quantities along the perturbed trajectory r'(t) are primed. The displacement 
of the element is defined as Sr = r' — tq. Keeping only linear terms in Sr we obtain the 
acceleration of the displacement, dr: 



Here, all pressure derivatives are evaluated along the unperturbed trajectory, (r = ro(t), t), 
and S is the entropy. 

In the limit of a non relativistic flow with a non relativistic equation of state Eq. (4) 
reduces to the equation derived using a somewhat different method by Goodman (1990). 
For an incompressible fluid, c s — >■ 00, in a steady state where pressure gradient supports the 
fluid against gravity, p~ 1 dp/dr = —g, Eq. (4) gives the classical Rayleigh- Taylor stability 
criterion, dp/dr < 0. 



Eq. (4) describes the behavior of the Lagrangian displacement of a perturbed fluid 
element. The growth of the displacement does not indicate, however, that a convective 
instability exists. Consider, for example, the flow of an expanding sphere where the velocity 
increases outward. If a fluid element is displaced from its original shell to some outer shell 
and continues to move with the outer shell, the displacement will grow although the fluid 
element does not overtake outer shells, i.e. although the instability is not convective. 
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3.3 A Qualitative Criterion for Local Convective Stability 
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Convective instabilities exist only if the acceleration experienced by the displaced fluid 
element is larger than the unperturbed acceleration of the fluid at the new position of the 
fluid element (rather than that at the original position of the element). The difference 
between the acceleration of a displaced fluid element and the unperturbed acceleration of 
the fluid in the displaced position at the rest frame of the shell into which the element had 
been displaced is 



da — 



c 2 dp 1 



e + p dr e + p 



de \ dp de 
dp J s dr dr 
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Sr, (5) 



(e + p) 2 \dr 

where the derivatives are evaluated at the shell rest frame and de/dp = [(de/ dr) / (dp/ dr)]. 
Using e = nmc 2 + aT and p = aT 4 /3, eq. (5) may be written in the form 
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where S, the entropy per baryon, is logT 3 /n. The subscript RF denotes that the deriva- 
tives are to be taken at the shell rest frame. 

Eq. (6) implies that convective instability will occur whenever the rest frame pressure 
and entropy gradients are parallel. This criterion is similar to the Schwarzschild criterion 
for convection is stars. However, at the case at hand, the fluid background flow is not 
stationary. Thus, in order for the instability to be of importance it is necessary that the 
time scale for instability growth be comparable to or shorter than the flow dynamical time 
scale. Only than a true dynamical instability will exist. We therefore compare the typical 



time in the observer frame for instability growth, 7~i ~ r y^/Sr/Sa 1 with the typical time for 
shell acceleration, the dynamical time r 7 , 
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The ratio of these times is 
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We conclude that: (i) A fireball flow is convectively unstable in regions where the en- 
tropy per baryon decreases outward (the pressure is assumed to decrease outward), (ii) As 
long as the fireball is highly dominated by radiation, rj 3> 1, the typical time for convective 
instability development is large compared to the dynamical time, 6>1. However, as the 
shell approaches the matter dominated phase, i.e. when rj ~ few, convective instabilities 
will develop rapidly. A rapid instability development may be obtained for larger rj values 
if the entropy scale length is considerably smaller than the pressure scale length. 

The displaced fluid elements are accelerated by radiation pressure. The energy driv- 
ing the instability is therefore the radiation energy. At the stage when the instability is 
expected to develop rapidly, i.e. when r\ ~ few, most of the shell energy is in the radiation 
field and had not been transformed to baryon kinetic energy. This energy is therefore 
available for driving the instability. In the limit of r] <C 1 most of the radiation energy had 
already been converted to baryon energy, and only a small fraction of it is available for 
the instability. At this stage the instability will continue to develop due to the energy it 
acquired during the transition phase from radiation to matter dominated evolution. How- 
ever, the energy associated with the instability will not increase substantially during this 
phase (The displacement will continue to grow, however without substantial acceleration). 

4. The Competition between Instability Growth and Shock Formation 

Shock waves are formed in the (Lagrangian) fireball regions where dr/i/dr < 0. Unless 
the energy density increases steeply with such regions will also have dS/dr < 0. This 
implies that regions where shocks will develop, i.e. where drji/dr < 0, are quite generally 
also convectively unstable. As explained in §2, a shock wave is formed at R sc = R v +i] 2 Ar 
as the inner high rj shells overtake the outer low t] shells (Ar is the initial separation of 
the shells). The instability, on the other hand, begins to operate at R ~ R v . Thus, if 

7 



Ar 3> Ri/rj then R sc — R v ^> R v and the shock forms long after the instability started 
developing. If, however, the inverse condition Ar<Ri/i] holds, instability development and 
shock formation are competing processes. The following numerical example demonstrates 
that even in this case the instability develops substantially prior to shock formation. 

We have calculated numerically the evolution of a fireball with an initial radius Ri = 
10 7 cm, initial temperature Tj = lMeV, and r\i profile decreasing outward from r\h = 50 to 
i]i = 5 over Ar ~ 7 • 10 5 cm. Fig. 1 presents the solution of (4) for the evolution of the 
displacement of a fluid element initially lying at the r\i transition layer (we have chosen 
5r > and 8r = at t = U). The displacement grows by the time of shock formation, 
r /Ri « 30, by two orders of magnitude. The instability is indeed convective, as seen 
by comparing the displacement evolution to the evolution of the distance between two 
adjoint shells (Note that the initial rapid displacement growth, obtained for r /-Ri<l, is 
not convective and is due to the differential expansion) . The results presented in fig. 1 are 
consistent with the qualitative analysis of §3.3, according to which fireball regions where 
dS/dr < are subject to substantial convective instability growth during the transition 
from radiation to matter dominated evolution. A shock wave forms only at the end of the 
computation shown in fig. 1. This demonstrates that the instability indeed won over the 
shock formation. It should be noted that as the inner shell overtakes the outer shell some 
of the baryon kinetic energy is converted back to radiation energy as the pressure in the 
shells interface region increases. This energy is then available for driving the instability. 

5. Conclusions and Applications to Astrophysical Fireballs and Grbs 

We have seen that a convective instability develops whenever dS/dr < 0. This in- 
cludes, quite generally, the important case of dr\jdr < 0, which appears whenever an inner 
high r\ shell accelerates an external low r] shell. The instability operates from the stage of 

8 



transition from radiation to matter dominated evolution, rj ~ few and R ~ R v , until the 
fireball becomes optically thin at R T (or at r^l 2 R T where the electrons decouple from the 
photons). The instability competes with shocks that form at shell crossing within AR SC 
after R v . 

We turn now to an astrophysical fireball which appears in a cosmo logical 7-ray burst. 
The energy release is 10 51 ergs and the rising time indicates that the initial radius Ri is 
around 10 7 cm (or less). With these parameters we have: R v = Rit]i = lC^cm/t^ry^ 
(where R i7 = ^/10 7 cm and rj i4 = WlO 4 ), Rr = (a T E /^rripC 2 ) 1 / 2 = 6 • 10 12 cm£?J( 2 7/^ 1/2 
(where E 51 = £/10 51 ergs) and AR SC « Arrf = 10 14 cmAr 6 7/ I 4 = 10 3 R v (Ar 6 /R i7 )rj H . 
Thus, for 77 < 3.3 • 10 4 El( 3 R~ 7 2 ^ s the instability will operate during a period when the 
fireball expands by a factor of ten or more. This gives ample time for the growth of the 
instability. Shell crossing and shocks occur within AR SC from R v . If Ar is not too small 
then AR SC ^> R v . However, even when this condition is not satisfied we have demonstrated 
that the instability operates rapidly enough before the shock forms. We conclude that 
generally, the instability has enough time to develop substantially whenever the conditions 
for its existence are fulfilled. 

We stress that our stability analysis presented is valid for a quasi-steady wind flow 
generated when the energy is released over a time scale larger than the light crossing 
time of the source, as well as for a fireball flow where the energy release is instantaneous. 
The (Lagrangian) regions of a wind created during periods when the emmited entropy 
increases with time, leading to a spatial profile with dS/dr < 0, will be unstable. This 
might have important implications to the recent model, suggested by Rees & Meszaros 
(1994), according to which GRBs are produced due to the formation of shocks between 
adjoint shells in the fireball. The formation of such shocks requires a wind where r\ increases 
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with time, resulting in regions with drj/dr < 0. Such regions are likely to satisfy dS/dr < 
(unless the source luminosity varies enormously with rj) and are, therefore, unstable. This 
implies that the effects of the instability must be taken into account within the framework 
of this model. 

The convective fireball instability described here will give rise to convection zones 
where high S bubbles penetrate into low S shells and vice versa. It is difficult to estimate, 
at present, the full implications of the instability to fireballs and to the observation of 
7-ray bursts, since for such an estimate it is necessary to consider the non-linear regime 
of instability development, where non-linear dissipation affects bubble motion. In the 
following we list several possible effects. First and most important is the fact that a 
situation where a high r\ shell accelerates a low 77 shell is unstable. The instability may 
therefore provide a mechanism for allowing high rj shells to penetrate through low 77 ones, 
thus allowing to achieve Lorentz factors higher than that corresponding to the average r\. 
This may ease the stringent limit on the allowed baryonic load. 

The mixing induced by the instability of high and low r\i shells will give rise to regions 
where the baryon density n fluctuates over length scales much shorter than both the 
fireball thickness and the source size, which are comparable and of the order of Ri. If 
the burst is produced by a collision with the ISM or by collisions of sub-shells within 
the fireball, this small scale structure may account for the rapid temporal variations of 
the observed pulse. This effect is especially important for the scenario where the burst is 
produced by a collision with the ISM. In this case, the duration over which an observer 
receives the radiation from a single sub-shell of the fireball is (Katz 1994) R c /rjiC, where 
the collision radius R c is (Piran 1994) R c = 1.3 • 10 15 cmK^ 3 'r]~ 2 ' 3 'nj 1 ' ' 3 (m is the ISM 
number density in units of cm -3 ). This duration is much larger than the source dynamical 
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time, R c /rjiRi = 1.3 • 10 4 i&[i 77 i4 . Thus, it would be hard to explain the rapid 

temporal variations of the pulse as resulting from temporal variations of the source. 

Finally, the generation of local random motions within the fireball shell will cause the 
fireball evolution to differ from that of a spherical smooth model due to the mixing induced 
by these motions and due to the kinetic energy associated with them. Such motions may 
further contribute to the build up of a magnetic field, which is necessary for reconverting 
the baryon kinetic energy to radiation energy. 
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Figure Captions 

Fig. 1: The evolution of a fireball sub-shell initially lying at the r\i transition layer. 
The horizontal axis measures the distance which the shell moved, which is proportional to 
the time that passes from the beginning. Solid line: The displacement 8r normalized to its 
initial value 5r;; Dashed line: ?y; Dash-dot line: 6; Dotted line: The distance between the 
shell and an outer adjoint shell normalized to the distance at t;. Note that a shock wave 
forms only at the end of this computation. 
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